Vacuum energy of two-dimensional M = (2,2) super Yang-Mills theory 
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We measure the vacuum energy of two-dimensional TV = (2, 2) super Yang-Mills theory using 
lattice simulation. The obtained vacuum energy density is fo = 0.09(9)(lg°)(;'^, where the first error 
is the systematic and the second is the statistical one, measured in the dimensionful gauge coupling g 
which governs the scale of the system. The result is consistent with unbroken supersymmetry, 
although we cannot exclude a possible very small non-zero vacuum energy. 
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I. INTRODUCTION 

Lattice simulation of gauge theory is of great im- 
portance since it provides a non-perturbative approach. 
For the supersymmetric systems, both theoretical for- 
mulation on the lattice and simulation techniques are 
now rapidly developing (see Ref. [l| for a review and 



Refs. ^,^,A,^,, 6, 7, 8,i9|,[lC 




for recent related 



works). One of the mysterious points of supersymmetry 
(SUSY) is that it is broken in our universe, although 
it is widely believed that SUSY is a relevant symmetry 
beyond the standard model. Because SUSY is such a 
strong symmetry, it is difficult to break it spontaneously 
in perturbative calculations, unless it is broken at the 
tree level. The lattice simulation will help to understand 
the non-perturbative effects for the SUSY breaking. 

In a recent work Q, on the basis of a formulation by 
Sugino the author together with Suzuki explicitly 
confirmed that lattice artifacts which break supersymme- 
try in fact disappear in the continuum limit. The target 
system was two-dimensional TV = (2, 2) super Yang-Mills 
theory (for other lattice formulations of this system, see 
Refs. HB, HI [13, [H, m, US [U ) . In order to suppress the 
fluctuation along the flat direction, we added scalar mass 
terms and numerically observed the Partially Conserved 
SuperCurrent (PCSC) relation, which is the conserva- 
tion law of the supercurrent modified by the scalar mass. 
Furthermore, we illustrated some physical applications, 
power-like behavior of certain correlation functions and 
static potential between fundamental probe charges [l3l| . 

The vacuum energy is the order parameter for SUSY 
breaking (2^ and we are interested in the spontaneous 
SUSY breaking of this system. Hori and Tong [2^ 
pointed out the possibility of spontaneous breaking while 
the previous simulation presented in Refs. [1, @ was con- 
sistent with no SUSY breaking. The vacuum energy ob- 
tained by the above simulation, however, had a rather 
large error because of an inefficient simulation algorithm. 

The purpose of this letter is to measure the vacuum 
energy with a small error which would help to determine 



whether supersymmetry is spontaneously broken or not 
in this system. We need, of course, to separate the spon- 
taneous supersymmetry breaking from a breaking caused 
by lattice artifacts. The result in Ref. shows that the 
latter does not survive in the continuum limit. 

The target two-dimensional system is a nice arena 
for simulation of supersymmetric systems, although it 
is a toy system from a viewpoint of phenomenolog- 
ically realistic four dimensional models. Since it is 
two-dimensional, computational cost is less than four- 
dimensional case. Furthermore, there is no sign problem 
in the continuum target theory, which also makes simu- 
lation easy. Important differences from four-dimensional 
TV = 1 supersymmetric Yang-Mills [28| are the existence 
of massless scalar fields and the existence of fiat direc- 
tions in the classical potential. 



II. FORMULATION 

We use the lattice model in Ref. and the method 
proposed in Refs [1, Q to measure the vacuum energy. 

The action of the target TV = (2, 2) super Yang- 
Mills theory is obtained by dimensional reduction from 
4-dimensional TV = 1 super Yang-Mills theory. In the 
Euclidean continuum spacetime, it is 



S=l I <PxiA \fmnFmn + ^^CTmDm^ 



(1) 
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Here, M and N run from to 3. Fm are the gamma ma- 
trices in four dimensions and C is the charge conjugation 
matrix. The covariant derivative Dm is Bm + *[^m,'], 
where 82 — 83 — 0, ^0 and Ai are two-dimensional gauge 
fields, A2 and are real scalar fields. The field strength 
is Fmn — 9m An — On Am + i[^A/,^Af]- The fermion 
^ = (V'Oi V'l? Xi '7/2)"^ is a 4 component spinor and H 
is a bosonic auxiliary field. The gauge coupling g has 
mass dimension 1 and we use it as a unit for dimension- 
ful quantities. The action is invariant under four super 
transformations and one can rewrite it in a supercharge 
exact form. 

The lattice action we use is ^lattice = -Ssugino + 
S'scaiar mass [3], whcrc S'sugino IS the lattice version of ^ 
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and has one exact nilpotent symmetry Q at finite lattice 
spacings out of four supercharges • "Sscaiar mass , which 
softly breaks the supersymmetry as well as Q-symmetry, 
is needed to lift the flat direction of the potential, be- 
cause otherwise the simulation docs not give thermalized 
configurations [25| . We thus define the system as an ex- 
trapolation in which the scalar mass fi goes to after 
taking continuum limit. In taking the continuum limit, 
we use fixed scalar mass in a physical unit. 

The vacuum energy is the order parameter for SUSY 
breaking. It is positive if and only if SUSY is sponta- 
neously broken. Therefore it is of crucial importance to 
choose the correct origin of the energy. Our choice is to 
use a Q-exact Hamiltonian (density). In Refs. 0, Q, the 
author together with his collaborators discussed that this 
choice is consistent with a property of the Wittcn index. 
They also confirmed that this method reproduces known 
results for supersymmetric quantum mechanics. For the 
current system, the Hamiltonian density is 7i = ^Q>Jq^\ 

where J^q^^ is the 0-th component of the Noether current 
corresponding to another supercharge Qq. Note that Q 
and Qq satisfy the super algebra {Q,Qo} — —2ido in 
the continuum. Since we do not have a lattice version 
of the Qq transformation, we use a naive discretization 
of the continuum current (see Ref. Q for its explicit ex- 
pression). We do not include the scalar mass term in the 
Hamiltonian density because the Noether current J''^"^ 
does not contain it. The Hamiltonian density is exactly 
identical to the one treated in Refs |4, 5]. 

We utilize the temperature as a conjugate external field 
to the Hamiltonian, i.e., we impose antiperiodic bound- 
ary condition (aPBC) for fermion fields in time direction. 
We use a finite spacial size Lg and temporal size /3, where 
f3 is the inverse temperature as well. What we measure 
is 

/oyx /aPBC d/i exp(-^iatticc) 
JaPBC exp(-S'iatticc) 

where we explicitly indicate the boundary condition, d/i 
denotes an integration measure of the path integral. If 
we adopted the periodic boundary condition, since the 
expectation value of Hamiltonian would be proportional 
to a derivative of the Witten index w.r.t the coupling 
constant, the expectation value would be always zero [2^ . 
Therefore, we could not obtain useful information for the 
SUSY breaking. We will take the zero temperature limit 
at the final stage to obtain the vacuum energy density: 
So = lim^^oo(linip^o(linia^o(H))) [13] ■ 

III. SIMULATION 

We evaluate the expectation value ^ using a Monte 
Carlo simulation. The detail of generating the configura- 
tions is found in Refs. 0, HI] [3l|. We use the Rational 
Hybrid Monte Carlo (RHMC) method together with the 
multi-time step for molecular dynamics. Compared with 
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FIG. 1: Distribution of the argument of the PfafRan for 200 
samples. The scalar mass is /i = 0.5g. The argument is cen- 
tered around 0, which justifies the replacement of the Pfaffian 
by the positive square root of the determinant. 



the method used in , which uses quenched configura- 
tions and rcwcighting method to treat the fermion effect, 
the statistical error is well under control. 

The effect of the complex phase of the Pfaffian of the 
Dirac operator is taken into account by phase reweight- 
ing of the square root of the determinant. Note that the 
target system in the continuum has real positive Pfaf- 
fian. Because our choice of the lattice spacing is small 
enough and the argument of the determinants is dis- 
tributed around zero, there is almost no sign ambiguities 
caused by taking the square root. In Fig. [T]and Fig. [21 
we plot the distribution of the argument of the determi- 
nant and a sample distribution from a direction calcula- 
tion of the Pfaffian ^] . FigJT] shows that the Pfaffian 
is almost real and positive [33| even for the largest lat- 
tice spacing and the lowest temperature, and guarantees 
the replacement of Pfaffian by the positive (i.e., the real 
part is positive) square root of the determinant. Note 
that the direct calculation of the PfafRan is numerically 
very expensive. From Fig. [2l one can see that as the lat- 
tice spacing becomes small the determinant and hence 
the Pfaffian approaches real and positive values, which is 
consistent with the property in the continuum. 

In the measurement, we calculate the Hamiltonian den- 
sity every 10 trajectories. We use binning technique (a 
typical bin size is 5) and use the jack-knife method to 
estimate the statistical error for each parameters. The 
continuum limit is taken as a linear extrapolation w.r.t. 
the lattice spacing a. The massless extrapolation is taken 
as a linear extrapolation w.r.t. /i^. We use fit and es- 
timate the statistical error of the fitting parameters as 
the value for which increases by 1. 

The common parameters are the following. We use 
SU(2) gauge group. Physical spacial size Lg is fixed to 
Lsg = 1.414. The parameter e for the admissibility con- 
dition in ^sugino is 2.6. 
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FIG. 2; Distributions of the argument of the determinant. 
The upper plot is one of the broadest distribution and the 
lower is one of the narrowest. For each plot, the scalar mass 
is = 0.5g. Plot with f3g — 4.243, ag — 0.236 uses the same 
parameter as the plot in Fig. [T] Even in the broadest cases, 
the argument is centered around zero. 



IV. RESULTS 
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FIG. 3: Example of the continuum limit, for /3g = 0.707 
and fi — 0.7g. We use data set with Pg = 0.707 to estimate 
systematic errors associated with the choice of the fitting func- 
tion and the fitting region. 
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FIG. 4: Examples of the extrapolation of the scalar mass fi 
to 0. The linear extrapolations in fT^ fit very well. 



We list the detailed parameters and raw results in Ta- 
ble HI For large /?, we use 3 different lattice spacings to 
obtain the continuum limit while for small /3 we use 4 
or more. It is because that in the former cases the ob- 
tained values are almost constant within the statistical 
error over the different lattice spacings; in the latter cases 
they depend on the lattice spacing. 

We estimate a systematic error in taking the contin- 
uum hmit using the data set with (3 = 0.707^ (Fig. [3]). 
Although we adopt the linear fitting to obtain the contin- 
uum limit, we tried a quadratic function A'a^ + B'a + C 
as well as a linear one Aa -I- i? in the lattice spacing a. 
After extrapolating /i — > 0, we obtained two different val- 
ues and we regard the difference as the systematic error 
in choosing the fitting function, which is 4.0%. To es- 
timate the errors associated with a choice of the fitting 
region, we repeated a similar analysis with linear fitting 
using small 4 lattice spacings and large 3 lattice spac- 
ings. Comparing these two results with the one obtained 
from all 6 lattice spacings, we estimate the error asso- 
ciated with the choice of the fitting region at 6.8%. In 
total, the systematic error in taking the continuum limit 
is 11%. 

After taking the continuum limit, we extrapolate the 



scalar mass ^ — > (Fig. |4]). As the figure indicates, a 
possible systematic error in the extrapolation is small 
compared with the statistical ones and thus is negligible. 

The final result is given in Fig. O We have fitted 
the result using two different functions: a power of /5, 
£{[3)/g'^ = ao{(3g)~'^^ + 02, and an exponential of /3, 
S{P)/g^ = 60 cxp(-6i/3g) + 62- We obtain gq = 1.941^^, 
ai = 2.59tl, ^2 = 0.09t^° and 60 = 157^}^, 61 = 
S.OltiJ^, 62 = 0.37tg. The errors are only statistical. The 
values of per degrees of freedom are 1.21 and 3.81 for 
7 degrees of freedom, respectively. Using the result of the 
fit by power function, we obtain — O.OOlg"^^ for the 
central value and the statistical errors. 

However, because the above fits use small (3 as well 
as large f3, one cannot exclude possible effects from the 
excited states. In fact, using only large (3 {(3g > 0.9), 
we obtain reliable fits for both exponential and power, 
which gives Sq = 0.17tn5^ and £0 = 0.05t.llg^, respec- 
tively [3J|. These values are consistent with the above 
central value within the statistical error. Using these 
values, we estimate that the systematic error in taking 
zero temperature limit is O.OSg^. 

What is the physical difference between the power and 
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FIG. 5: Energy density £{l3). Errors bars denote statistical 
errors only. The dashed curve is the best fit of data with a 
power function of (3. As the inverse temperature (3 becomes 
large, the energy density approaches 0. 

the exponential fitting? If the spectrum has a gap, the 
energy behaves as an exponential function of /3. A power 
behavior of (3 comes from gap- less excitations [s^. In a 
finite volume system, usually we have a discrete spectrum 
which leads to a mass gap. However, the current system 
has the flat direction in the potential at least in the clas- 
sical level, and the fluctuation along the flat direction 
can provide a continuum spectrum. The current result 
does not allow us to determine which is the case [s^. If 
the flt by the power was preferred, it would be an evi- 
dence of the recovery of the flat direction after the /i ^ 
extrapolation. 

In total, we obtain the vacuum energy £q = 



0.09(9) (^g");?^, where the first error is the systematic and 
the second is the statistical one. It is consistent with 
within the error, that is, it is consistent with supersym- 
metry not being spontaneously broken. 

V. CONCLUSION 

In this letter, we measured the vacuum energy of two- 
dimensional TV = (2, 2) super Yang-Mills theory using 
lattice simulation. Because our simulation does not have 
explicit supersymmetry breakings caused by lattice arti- 
facts, the vacuum energy we measured works as an or- 
der parameter of spontaneous SUSY breaking. The re- 
sult is consistent with within the error, which seems to 
imply that supersymmetry is not spontaneously broken, 
although we cannot exclude the possibility of a non-zero 
value smaller than the error which would indicate a spon- 
taneous breaking of SUSY. 
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TABLE I: Expectation value of the Hamiltonian density 1(H) and number of the independent configurations after binning for 
each parameters, g is the dimensionful gauge coupling. 
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